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Abstract 

We present new predictions for the total cross section of squark pair-production at Tevatron and 
LHC through next-to-next-to-leading order within the Minimal Super symmetric Standard Model. 
The results are based on the numerically dominant soft corrections. They are exact in all logarith- 
mically enhanced terms near threshold, include the Coulomb corrections at two loops and exact 
scale dependence. We translate the increased total cross section at next-to-next-to-leading order 
into improved exclusion limits for squark masses and we investigate the scale dependence as well 
as the sensitivity on the parton luminosity. 



1 Introduction 



Supersymmetry offers an attractive way for possible extensions of the Standard Model. Its most 
popular incarnation for phenomenological studies is the Minimal Supersymmetric Standard Model 
(MSSM) which features a rather rich spectrum of new (heavy) particles. Active searches are cur- 
rently being performed at the Tevatron and soon at LHC. Squarks and gluinos as scalar and spin- 
1/2 supersymmetric partners of quarks and gluons could be pair- produced at a hadron collider at 
sizable rates (assuming 7?-parity). So far, the Tevatron (CDF and DO collaboration [1-3]) has been 
providing lower limits on their masses depending on certain assumptions about the parameter space 
of the MSSM. At LHC in contrast, one expects for an initial luminosity of about lOfb 1 and typi- 
cal values of the MSSM parameters 10 6 (2500) squark pairs for a mass of 300 GeV (1 TeV) [4,5]. 
Given these large rates and the dedicated searches performed, it is natural to investigate the accu- 
racy of the available theoretical predictions including quantum corrections. 

With squarks (and gluinos) carrying color charge, it is not surprising that Quantum Chromo- 
dynamics (QCD) provides the dominant corrections to the production cross section. This fact was 
realized some time ago and led to the computation of the complete next-to-leading order (NLO) 
QCD corrections [6]. The upshot is a large increase of the rate in comparison to leading order 
(LO) QCD predictions along with a reduced scale dependence indicating the improved theoretical 
uncertainty. The origin of these large higher order QCD corrections in hadro-production of (heavy) 
colored particles is well known, since it is related to universal QCD dynamics. As a typical pat- 
tern so-called Sudakov logarithms show up which originate from soft gluon emission in regions of 
phase space near the partonic threshold. They depend on the squark velocity P = (1 —4m?/,?) 1 / 2 

and become large for center-of-mass energies near the threshold for squark pair-production, 
s ~ 4m~. Sudakov logarithms can be organized to all orders of perturbation theory by means of a 
threshold resummed cross section to a given logarithmic accuracy. 

In this letter, we study soft gluon effects for the total cross section of squark hadro-production. 
We employ Sudakov resummation to generate approximate next-to-next-to-leading order (NNLO) 
QCD predictions which are accurate in all log(P)-enhanced terms at two loops. Moreover, we 
include the complete two-loop Coulomb corrections as well as the exact dependence on the renor- 
malization and factorization scale. To that end, we largely follow a similar study for top-quark 
hadro-production at the Tevatron and the LHC [7, 8]. The importance of Sudakov resummation for 
the latter reaction has often been emphasized in the literature (see e.g. Ref. [9]). Squarks are gener- 
ally believed to be heavier than top quarks, but light enough to be pair-produced at these colliders 
as well, so that the need for threshold resummation (and the technique) carries over from top-quark 
pair-production. Recently, the soft gluon resummation for squark and gluino hadro-production has 
been performed to next-to-leading logarithmic (NLL) accuracy in Ref. [10], and results compati- 
ble with ours were found. See also Refs. [11-13] for the electroweak contributions to (top)-squark 
pair-production through NLO. 

The letter is organized as follows. We recall the dominant parton channels contributing to the 
cross section of squark pair-production. Then we discuss the NLO QCD corrections. For the latter, 
we have determined fits to the exact NLO scaling functions for representative choices of squark 
and gluino masses (presented in the Appendix). Next, we describe the steps necessary to achieve 
soft gluon resummation to next-to-next-to-leading logarithmic (NNLL) accuracy where, we focus 
on the differences with respect to the well-known procedure for top-quarks. This extends the NLL 
results of Ref. [10]. We employ the resummed cross section to derive new (approximate) NNLO 
expressions for the scaling functions and analyze the dependence of the hadronic cross section on 
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the scale and on parton distribution functions (PDFs). For the finite order expansion to NNLO, 
we find good apparent convergence properties and a markedly improved stability of the total cross 
section with respect to scale variations. We give two examples, how our results translate into new 
limits on the squark masses. 

2 Setting the stage 

We focus on the inclusive hadronic cross section of hadro-production of squark pairs, o pp ^qq*x 
which is a function of the hadronic center-of-mass energy y/s, the squark mass niq, and the gluino 
mass m.g. In the standard factorization approach of perturbative QCD, it reads 

a pp ^ W x(s,ml,mj) = ^ J ds Lij(s,s,/jj) d i7 ^*(£,m?m? , (1) 
where the parton luminosities L ;; - are given as convolutions of the PDFs fi/ p defined through 

s 

Lij(s,s,nj) = ±J^f i/p ^}^y j/p (^ 2 f ,-J, (2) 

s 

where s denotes the partonic center of mass energy. Factorization and renormalization scales [Xj 
and /j f . are identified (jjr = /j r = //). The sum in Eq. (OQ) runs over all massless parton flavors 
(«/ = 5 for LHC and Tevatron). The hard parton scattering cross section dij^qq* appearing in 
Eq. (OQ) receives at Born level contributions from the channels 

88 -> qktii, (3) 

qiq~i -> qkqh ( 4 ) 

where /, j, k and / denote flavor indices. For the calculation of the hadronic cross section (OQ) in this 
paper, we always sum over all final state squark flavors allowed by quantum number conservation 
in proton-(anti)-proton scattering except for top-squarks (stops). The contributing diagrams are 
displayed in Fig.Q] 

While the gluon fusion process © always leads to like flavor squark pairs (k = I) and identical 
chiralities (see e.g. Ref. [6] for details on the quantum numbers), a qq initial state leads to a 
more interesting flavor structure of the final state squark pair. For the qq- scattering process © 
we have final states with the same flavor structure if the flavors of the initial state are different 
(i^j^i = kAj = I) due to a gluino exchange in the ^-channel. Likewise, if the initial flavors 
are equal (i = j), then the final flavors may be different from the initial flavors, but the flavors 
of the qq* -pair are equal due to a gluon exchange in the s-channel (i = j => k = I). Moreover, 
if we keep the squark chirality as a second quantum number, we encounter a richer structure and 
much of the discussion (including soft corrections) carries over. For the subsequent study of NLO 
and approximate NNLO perturbative QCD, we restrict ourselves to cross sections summed over 
all possible final squark flavors and chiralities (L,R). The complete list of Born cross sections for 
different flavors and chiralities for the process (HI) is presented in the Appendix. In the following, 
we focus on the total partonic Born cross sections 6 s where all flavors and chiralities are summed 
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Figure 1: Contributing diagrams to the processes gg — > qq* (diagrams (a)-(d)) and qq — ► qq* (dia- 
grams (e),(f))- Diagram (f) is proportional to the squark-gluino Yukawa coupling a s and introduces de- 
pendence on the gluino mass nig. 



over, using the conditions i = j or i ^ j (see e.g. Ref. [6]), 
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with the strong coupling constant oc 5 and the squark-gluino Yukawa coupling & s identical to a s as 
required by supersymmetry (see the discussion below). Moreover, we have abbreviated 



log 



s + 2(m~ — mf) -5(3 
f+2(m|-m?)+fp 



(7) 



where (3 = (1 — 4m~/s) */ 2 is the squark velocity. 

For the studies of higher order QCD corrections, the partonic cross section can be expressed in 
terms of scaling functions fy. For gluon fusion, we define 



Ggg^qq*{S,>nq) 



(«), 



J f V 



1=0 



(8) 



and, likewise, for gg-scattering fj£9 with additional dependence on the gluino mass nig, see Eq. ©. 
Here T| = s/ (4m?) — 1 is a measure for the distance to the production threshold at s = 4m? The 
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NLO QCD corrections are known [6] although they are not available in analytical form (unlike 
the case of top-quark hadro-production [14]). Instead in Ref. [15], the authors present a numerical 
program called Prospino to calculate the complete hadronic cross section at NLO. We have used 
this program to extract values for the scaling functions fgl°\ fjfi°\ and (all terms proportional 



to log (ju 2 /m~) will be discussed below). Subsequently, we have determined fits of these functions 
to per mille accuracy based on the following ansatz, 
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The threshold logarithms log(P) as well as the Coulomb corrections (~ 1/(3) are kept exactly [6], 
while the fit function /i(p) is given in Eq. (IA.12I) . Near threshold, /i(P) = a\$ + ((3 2 ) holds (see 
Tab. [3]) and these fits allow for an easy handling of the NLO QCD corrections in phenomenological 
applications. 

The log(P) terms appearing in Eqs. © and (flOl) can be resummed systematically to all orders 
in perturbation theory employing the well established techniques (see e.g. Refs. [16-18]). The 
resummation typically proceeds in Mellin space after introducing moments N with respect to the 
variable p = 4m~/s as 



o{N,m%m\) 



&(s,mg,m Sj 



(11) 



For scattering reactions with non-trivial color exchange, one has to choose a suitable color basis 
for the total cross-section. It is convenient to select a decomposition according to color-singlet and 
color-octet final states, 
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and the octet terms d^- g can be easily derived from Eqs. © and ©. The resummed cross sections 
(defined in the MS-scheme), for the individual color structures of the scattering process are then 
obtained as single exponentials in Mellin- space, 



&ij,i(N,mjj,mj) 
dg r (iV,m2.,m|) 



8u,M 2 g,m 2 g )-exp {G ij:I (N+ 1)) + (AT 1 log' 7 N) , (15) 
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where all dependence on the renormalization and factorization scale n r and fir is suppressed and 
the respective Born term is denoted afj r The exponent G ;; .i contains all large Sudakov logarithms 

log k N and the resummed cross section (fT51) is accurate up to terms which vanish as a power for 
large Mellin-/V. To NNLL accuracy, G^i is commonly written as 

Gij^N) = logN ■ g}j(X) +g} JiI {k) + ^gliW +■■■, (16) 

where X = Po logiVoCj/ (47l) . The functions gy for the singlet and octet color structures are explic- 
itly given in Ref. [7] and can be taken over from the case of top-quark hadro-production (see also 
Ref. [10] for the results to NLL accuracy). All gf,, k = 0, ...,3 depend on a number of anomalous 
dimensions, i.e. the well-known cusp anomalous dimension A, the functions D and D qq control- 
ling soft emission, and the coefficients of the QCD [3-function. At higher orders their precise 
expressions also depend on the chosen renormalization scheme, thus on the dynamical degrees 
of freedom. At the center-of-mass energies of Tevatron and LHC and for the mass ranges cur- 
rently considered in MSSM phenomenology m q ,mg ~ 0(2OOGeV — ITeV), a scheme is appro- 
priate which decouples all heavy particles (top-quark, squarks, gluino). Thus, we are left with 
the Standard Model [3-function coefficients |3o = 11 — (2/3)rif and Pi = 102 — (38/3)ny and the 
same expressions for the anomalous dimension A, D and D qq as in the case of top-quark hadro- 
production (see Ref. [7]). Our scheme choice is in line with the exact NLO QCD calculation [6] to 
facilitate matching (see below). Throughout the paper, m q and mg denote pole masses. 

Having Eq. (TT3T ) and all quantities necessary for its explicit evaluation at our disposal, we use 
the resummed cross section as a generating functional for the threshold approximation to the yet 
unknown NNLO QCD corrections, i.e. the inclusive partonic scaling functions j^^ in Eq. ©. 

Substituting all numerical values and setting n/ = 5, we obtain in the MS-scheme the approximate 
NNLO results, 
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where the coefficients of the log 4 |3-terms are exact and the constants a\ from the NLO matching 
are given in Tab.[3l respectively Eq. (IA.14I) . In the derivation we have used the linearization of the 
Born functions, 



192 
4 

'mj 9^ (m 2 + m 2 )2 



fS = \t +0(V 3 ). (20) 



For our phenomenological studies though, we have always substituted the full Born result for f-j ^ 
in Eqs. dUb-dHb. 

A few comments on Eqs. (fT7l)-(fT8l) are in order here. The results are accurate to all powers 
in log(|3) at two loops. This has been achieved by consistent matching of the resummed cross 
section to the exact NLO result. To that end, we have used our fits © and (flOl) . specifically the 
constants a\ of Eq. (IA.12I) as given in Tab. [3l These constants enter in the coefficients of log 2 ((3) 
and log(P) in Eqs. (fT7l)-(fT8~l). i.e. the quadratic and linear logarithm, so that the first two terms 
become ( — 1309.5423, 3523.9200) in Eq. (TTTI) . Likewise for a representative choice of squark and 
gluino masses mq = 400 GeV and nig = 500GeV, we find the first two terms of log 2 (P) and log(P) 
to combine into (567.40606,932.92034), if the quarks in the initial state have identical flavors, 
and into (450.41 183, 923.62614) in Eq. (fT8l) , if the quarks in the initial state have different flavors, 
see also Eqs. ([201) and (1A. 14b . Strictly speaking, for the NLO matching separate constants are 
required for the singlet and octet color structures at NLO and for this reason the exact numerical 
coefficient of the term linear in log |3 will differ slightly. However, experience from studies for top- 
quark hadro-production shows this effect to be marginal [14] and well covered within our quoted 
overall uncertainty for the threshold approximation. On top of the threshold logarithms we have 
also added in Eqs. (fT7l)-(fT8l) the complete two-loop Coulomb corrections (as summarized e.g. in 
Ref. [7]). Finally, we comment on the presence of the squark-gluino Yukawa coupling a s , which is 
identified with the strong coupling constant a s as required by supersymmetry. In the MS-scheme 
this is achieved by a (finite) renormalization, which is correctly implemented at one-loop level by 
our matching procedure at NLO. At the two-loop level, the additional renormalization of a s would 

fn) (J) 

only affect the constant terms C q q and C q .q. in Eq. (PTSl . Since we have no control over these 

(2) 

constants, we set all two-loop constants to zero in our phenomenological studies. 

At this stage, it only remains to discuss those terms in Eq. ([8]) which describe scale dependence, 
i.e. the terms proportional to log(/j 2 /m~). Through NNLO, it concerns f^qj, f^qj and f^- 2 } which 
are entirely determined by renormalization group arguments (see. e.g. Refs. [19,20]). They can 
be constructed with the help of lower order results, that is the scaling functions and and 
the splitting functions Pjj. The latter quantities govern the PDF evolution. They can be expanded 
as 



(k) 

and explicit expressions for the P-j can be found in Refs. [21,22]. Following Refs. [19, 20], 
we can calculate the scale dependence of the partonic cross section. Again, we have identified 
renormalization and factorization scale n r = /jf = /j. However, the factorization scale dependence 
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can easily be recovered by re-expansion of OL s (/jf) in terms of CL s (j* r ). This yields the following 
result 



fij ,) = j^^r-C^-fP®^)' (2» 
(2i) _ 1 /or. ,(oo) .(oo)_„(i) ,(oo)^ D (i; 



kj 



^MP-C^-f^^)- C23) 

+3^ - |(WT - | Wf> ®n ( , 0) ) • <24) 

where <g) denotes the standard Mellin convolution, i.e. ordinary products in Mellin space under 
transformation CCD). The summation over repeated indices for admissible parton contributions is 
implied, although for phenomenological applications we restrict ourselves in Eqs. (1231) and (124)) to 
the (numerically dominant) diagonal parton channels at two loops. Note, that the scale dependent 
two-loop contributions are exact at all energies also away from threshold. In this context, we 
remark that the usual definition of the splitting function P qg (see e.g. Refs. [21,22]) needs to 
be divided with a factor 2rif, because the standard definition implicitly contains a sum over all 
flavors. One other comment concerns the two-loop splitting function P qq in Eq. (123T ). Obviously, 
only the flavor non- singlet splitting function P qq = P+ enters the convolution ([231 for the flavor 

(21) 

non-singlet scaling function f„ = because flavor number is conserved. Likewise for the flavor 
singlet case f^j , the flavor singlet splitting function P qq = +P ps (i.e. the sum of non-singlet 
and pure-singlet) contributes. Finally, f^j depends only on the ratio of the gluino and squark 

(\\\ (22) 

mass r q = nig/m q . Therefore through Eqs. (1221) and (1241) . the scaling functions f qjq -. and f q . q . 
depend only on this ratio. 



3 Numerical Results 

Let us start to illustrate the phenomenological consequences. In our numerical illustration we 
assume throughout this paper for the gluino mass the relation nig = 1 .25m q . Here our choice of pa- 
rameters m q = 400 GeV and nig = 500 GeV as a reference point has been influenced by the current 
limits from direct searches. Thus far, lower limits on m q and nig are provided by Tevatron (CDF 
and DO collaboration [1-3]) and an absolute lower limit of 379 GeV and 308 GeV, respectively, for 
the masses of squarks and the gluino in the mSUGRA framework (Aq = 0,/j < 0,tan[3 = 3) has 
been quoted. Moreover, the non- stop- squarks are expected to be in a narrow mass range, because 
the Standard Model quarks are (nearly) massless and, in this case, nearly mass -degenerate squarks 
are a property of mSUGRA. In the mSUGRA framework, the t\ is lighter and the ?2 heavier than 
the other squarks and a typical particle spectrum can be found in Ref. [23]. As announced already 
above, we have excluded stop quarks from our considerations here. Currently, the NLO QCD 
corrections to stop production are known [24] (see e.g. Ref. [25] for a search for the lightest stop 
fi). 
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Figure 2: Scaling functions (left column) and /^v (right column) with i = 0. . .2, j < /, for m f j = 
400GeV and m« = 500GeV. 
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We display in Figs. [2] and [3] our results for the scaling functions as defined in Eq. © up to 
second order in a s and separated according to the parton initial state gg and qq. The results for 



J 'j 



(20) 



(21) (2?) 

fij and ffj axe new. To illustrate the effect of threshold logarithms, we compare the scale 



independent functions fjj , 



and f;j directly in Figs. |2(e)[ |2(f)j and |3(c)j One clearly sees 



the effect of the large logarithms in (3 at low r\ giving rise to large perturbative corrections. At 
larger r\, the results for from Eqs. (fT7l)-(fT8l) vanish quickly and our approximation ceases to 
be valid for r\ > o(1...10). 
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Figure 4: The total hadronic cross section for pp — ► qq* at LHC with yfs = 14 TeV to LO, NLO, and NNLO 
in QCD as function of the squark mass mq. For the gluino mass we assume the relation mg = 1.25m f7 . 
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Figure 5: Same as Fig.@]for Tevatron with *Js = 1.96 TeV. 

Let us next investigate the consequences for the hadronic cross section CD for squark pair- 
production at LHC and Tevatron. The necessary convolution of the parton scaling functions with 
the parton luminosity © emphasizes the threshold region of phase space. The parton luminosity 
Ljj is steeply falling with increasing energies. Consequently, the total hadronic cross section is 
effectively saturated from partonic processes close to threshold with the kinematics being very 
similar to the case of top-quark hadro-production [7]. As an upshot, one can conclude, that our 
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approximate NNLO result captures the numerically dominant part of the complete (yet unknown) 
NNLO corrections. Eqs. (fT7l)-(fT8l) should thus represent a very reliable estimate. 

In Fig.|H we present the total cross section at the LHC as a function of the squark mass niq. We 
use the PDF set CTEQ6.6 [26], if not stated otherwise. In the left figure, we show the LO and NLO 
QCD predictions for the cross sections along with their error bands due to the scale uncertainty for 
scale choices x = /u/niq — 2 and x — 1/2. In Fig. @]on the right, we show the same comparison for 
the LO and NNLO prediction. While the error band of the LO cross section is quite large (due to 
the scale uncertainty of the strong coupling constant a s and the PDFs) it shrinks when considering 
the NLO QCD prediction. For the cross section at NNLO accuracy, we apply our newly derived 

(20) (21) (22) 

scaling functions f-j , jy and f-j of the previous Section and observe how the error band 
contracts to a thin line. This proves the significant reduction in scale uncertainty at this order of 
perturbation theory. Numbers for the hadronic cross section for different squark masses and scales 
are listed in Tab. [Q The NLO cross section increases by about 50% compared to the LO cross 
section, and likewise, the approximate NNLO cross section by about 9% compared to the NLO 
result. The latter value, i.e. 9%, is to be contrasted with the rather small ^-factor for the NLL 
resummed cross section reported in Ref. [10] of about 1 — 2% in the range niq = 200-1000GeV. 
As a direct phenomenological consequence we can translate our result into shifts of the exclusion 
limits for the squark masses for one example. For instance, a cross section of 30pb corresponds at 
LO, NLO, and NNLO to a squark mass of about 397 GeV, 430 GeV, and 437 GeV, respectively. 

In Fig. |51 we present the same plot for the total cross section of pp — > qq* at the Tevatron. 
Here, the cross section is more than two orders of magnitude smaller than at LHC due to the 
lower center-of-mass energy of Tevatron. Again, we observe improved stability of the perturbative 
predictions with respect to the scale variation. In Tab. |2l we give explicit values for the cross 
section and the NNLO result implies new exclusion limits. For example for a cross section of 
60 fb one comes up with a squark mass of about 370 GeV, 391 GeV and 396 GeV at LO, NLO, and 
NNLO, respectively. 

In Fig. [61 we show the contributions of the different parton channels © and © to the total 
cross section. At the LHC, the contribution of the gluon channel is dominant only up to a squark 
mass of about 370GeV. The ratio of the two channels, o(^2qq — > qq*)/o(gg — > qq*), changes 
from 45%/55% for m q = 300 GeV to 67%/33% for m 9 = 800 GeV. This is not what one naively 
expects from a p/?-collider. The reason for this behavior is the summation over all possible initial 
quark flavors and final squark flavors and chiralities. The gluon channel has only 10 contributing 
subprocesses, the quark channel however 140 subprocesses. At the Tevatron, the gluon channel is 
nearly negligible as one can see from Fig. |6(b)| and the large number of contributing subprocesses 
in the quark channel is additionally enhanced by the large ^-luminosity. The ratio of the two 
channels increases from 93%/7% for niq = 200 GeV to 98.5%/1.5% for m 9 = 500GeV. 

In order to illustrate the improved theoretical prediction due to the reduced scale dependence, 
we show in Fig. |7(a)| the LO, NLO, and NNLO predictions for the cross section at LHC with 
y/s =14 TeV. One observes a decreasing dependence on /j with increasing order of perturbation 
theory. Within the commonly chosen scale interval x = 0.5 . . .2, the LO, NLO, and NLO cross 
section varies between 37 pb and 23 pb, 45 pb and 40 pb, and 44 pb and 47 pb, respectively (see 
Tab. [T]). This amounts to a residual theoretical uncertainty of only ±3% for our approximate 
NNLO prediction and improves the value for the resummed cross section to NLL accuracy of 
about ±10% quoted in Ref. [10]. 

The presently discussed schedule for the initial phase of LHC includes operation at a center- 
of-mass energy lower than 14TeV, a value of lOTeV has often been quoted. For this reason, it is 
interesting to study the energy dependence of the cross section. In Fig. |7(b)| we display our new 
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Figure 6: The contributions of the channels gg — ► qq* and ^qq ^ qq* t0 tne tota l hadronic cross section 
at NNLO for the LHC (left figure) and the Tevatron (right figure). The gluino mass is trig = l.25rriq and the 
scale p = rriq. 
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Figure 7: The total hadronic cross section for pp — > qq* at LHC as a function of the scale p. (left figure) and 
as a function of the center-of-mass energy ^fs (right figure) for the reference point. On the left hand side we 
compare the LO, NLO, and NNLO predictions for the p dependence with y/s = \A TeV. On the right we 
display the ^/s dependence of our NNLO prediction at scale p = rriq. 



approximate NNLO result as a function of the center-of-mass energy. We observe, for instance, at 
10 TeV a cross section of about 19pb for our reference point {rriq = 400 GeV and rrig = 500 GeV) 
compared to 47 pb for the design energy of 14 TeV. 

The final point of interest is the sensitivity to the parton luminosity and the associated uncer- 
tainty. In Fig. [81 we show the NNLO cross section for squark pair-production together with the 
error band due to the PDF uncertainties. We used two sets of PDFs: CTEQ6.6 (left figure, [26]) 
and MRST2006nnlo (right figure, [27]). The PDF uncertainty has been calculated accordingly to 
Ref. [26]. At higher squark masses, the error on the cross section is significantly increased be- 
cause the relevant phase space probes the gluon luminosity in the high-x region where the gluon 
PDF has a large uncertainty. The relative error increases from ±3.5% at m q — 300 GeV to ±9% at 
rriq = 800 GeV. Note that the quoted error for the MRST2006nnlo PDF set is significantly smaller 



12 




m- [GeV] m- [GeV] 



Figure 8: The PDF uncertainty of the NNLO cross section for the two PDF sets CTEQ6.6 [26] (left figure) 
and MRST2006nnlo [27] (right figure). 



[GeV] 


o 


(LO)[pb] 
x=l 


x = 2 


o(NLO)[pb] 


x = 2 


o(NNLO)[pb] 

JC — ^ — 1 — 2 


300 


148.4 


116.6 


93.2 


179.5 


174.2 


161.9 


177.4 


189.8 


193.4 


400 


37.1 


28.9 


23.0 


44.6 


43.1 


39.8 


44.1 


46.8 


47.4 


500 


11.9 


9.3 


7.3 


14.4 


13.8 


12.7 


14.2 


15.0 


15.1 


600 


4.51 


3.49 


2.75 


5.47 


5.22 


4.78 


5.39 


5.66 


5.68 


700 


1.91 


1.47 


1.16 


2.33 


2.21 


2.01 


2.29 


2.40 


2.40 


800 


0.87 


0.67 


0.53 


1.07 


1.01 


0.92 


1.06 


1.10 


1.10 


900 


0.43 


0.33 


0.25 


0.53 


0.50 


0.45 


0.52 


0.54 


0.53 


1000 


0.22 


0.17 


0.13 


0.27 


0.25 


0.23 


0.27 


0.28 


0.27 



Table 1: Numerical values for the squark pair-production cross section at LHC with y/s = 14TeV and the 
CTEQ6.6 PDF set [26]. The QCD predictions are given at LO, NLO, and NNLO accuracy and for different 
squark masses and scales x = fi/rriq. 



than for the CTEQ6.6 set. Its error increases from ±1.5% to ±4% for squark masses ranging form 
300 GeV to 800 GeV. 



4 Conclusion and Summary 

In this letter we have investigated the effect of higher order soft corrections on the total cross 
section for hadronic squark pair-production. These radiative corrections make up numerically 
for a large part of the higher order QCD effects. Starting from the existing NLO calculation we 
provide for the ease-of-use the NLO scaling functions in the form of parameterizations which 
are accurate at the per mille level. Subsequently we have employed well established techniques 
for soft gluon resummation to derive new approximate NNLO results. Our two-loop expressions 
for the scaling functions fjj are exact in all logarithmically enhanced terms near threshold and 

(21) (22) 

they include the Coulomb corrections. All two-loop scaling functions /■ and /■ governing 
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(NNLO)[fb] 


[GeV] 




x=l 


x = 2 




x = 1 


x = 2 




x= 1 


x = 2 


200 


10735.0 


7641.1 


5606.6 


13394.7 


11836.7 


10091.7 


12707.5 


12749.7 


12100.4 


250 


2450.5 


1726.6 


1254.3 


3182.5 


2766.6 


2328.8 


3043.4 


3019.7 


2835.0 


300 


602.4 


419.9 


301.9 


816.9 


698.6 


580.4 


789.9 


775.2 


720.1 


350 


151.8 


104.6 


74.4 


215.2 


181.1 


148.4 


211.0 


204.9 


188.2 


400 


37.9 


25.8 


18.2 


56.3 


46.6 


37.7 


56.2 


53.9 


48.9 


450 


9.2 


6.2 


4.3 


14.3 


11.6 


9.2 


14.5 


13.7 


12.3 


500 


2.1 


1.4 


1.0 


3.4 


2.7 


2.1 


3.6 


3.3 


2.9 



Table 2: Same as Tab. Q] for Tevatron at y/s = 1.96TeV. 



the scale dependence have been computed exactly using renormalization group arguments. 

For Tevatron and LHC, our approximate NNLO cross section leads to a cross section increase 
of 9% compared to the NLO predictions, which translates into higher exclusion limits for squark 
masses. Moreover, with our approximate NNLO result we have found significantly improved 
stability with respect to variation of the renormalization and factorization scale (keeping /j r = /jy)- 
This leads to a residual theoretical uncertainty of 3% plus a (largely uncorrelated) error due to the 
parton luminosity depending on the particular PDF set. 

As a possible extension as far as the study of QCD corrections is concerned we would like to 
mention bound-state effects for squark pair-production at hadron colliders, which depend on the 
particular color and angular momentum quantum numbers of the qq*-pak. This would allow the 
resummation of the Coulomb corrections to all orders. From similar recent work for top-quark 
pairs [28,29] we would expect a shift of the total cross section by O (1%) due to these bound state 
corrections. 
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Appendix A: Useful formulae 

For the partonic cross section one can distinguish the following subprocesses for squark pair- 
production according to the initial flavors and chiralities (see also Refs. [13,32,33]): 



qq -> q R q* R , q =u,d,s,c,b, (A.l) 

qq — ¥ <?^^1' 9 =u,d,s,c,b, (A.2) 

qq QrQr, q,Q = u,d,s,c,b, q=£Q, (A.3) 

qQ -> qRQ* R , q,Q = u,d,s,c,b, q^Q, (A.4) 

qQ qR&t, q,Q = u,d,s,c,b, q^Q. (A.5) 
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The same processes are possible if one replaces the chirality index R with L and vice versa. The 
process (IA.1I) proceeds via gluon exchange in the ^-channel and gluino exchange in the /-channel, 
the process (IA.3I) only via gluon exchange, and the processes (IA.2I) . (IA.4I) . and (IA.5I) only via 
gluino exchange. The corresponding Born cross sections are given by 



G(qq -> <?^) 



27f P 



27ra 2 
21s 



6(3 + 3 



2(mj-m% ) 

1+ { 8 , qR> \u 



+ 



2%a s dL s 
21s 



P 1 + 



+ 2- 



2 \2 i 2~ 



, (A.6) 



27iaf 



2 ml A/K s + ^i-^) 2 -4m? ff m2 



9L 



qR qh 



9 s (mj-mj L )(mj-mj p )+mp 



qn' 



3 

275 P 



2^af 



(A.7) 



(A.8) 



9£ 



(m i +m L~ f)2_4m i m L 



+ i+ 



2m? 



L lR 



m 



Or 



o{qQ^q R Ql) 



(A.9) 



(A.10) 



9 s (ml-m^)(mj-ml g )+mj-s 
with m.q, m,q R/L , and »tg denoting squark masses, L\ being defined in Eq. (Q, and L\ given by 



4 



log 



f+2m|-m| fi -m? s 



m? + m? — f) 2 — 4m? m? 

1r Q R > C 1R Q R 



s + 2m? - m| s - m ? s + J {m\ R + m\ R - s] 



2 — 4m? m? 

9* Qr, 



(A.ll) 



We have parameterized the numerical result of Prospino [15] for the NLO scaling functions 
in fgl°\ fq q ; , and by a fit (accurate at the per mille level), where the fit function h($) in 
Eqs. © and (Oil reads 

fc(P) = pfai+a 2 p + a3p 2 log(8p 2 )+fl 4 (3 3 log(8p 2 ) 



1 



:log 2 (l+ri)+fl 6 - 



1 



:log(l+ri)+a7-— J— log 



x/T+rf ov v^+ri" DV ' ' ' "' 1+ri 

+a 8 zlog 2 (ri) +a 9 zlog(ri) + ai z 2 log 2 (r|) +anz 2 log(r|) 
+ai2Z 3 log 2 (ri) +ai 3 z 3 log(ri) +ai 4 z 4 log 2 (ri) +ai 5 z 4 log(ri) 

+ai 6 z 5 log 2 (r|) +ai 7 z 5 log(ri) 



1 + P 
1-P 



(A.12) 
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and z is defined as 



z = 



(i+ti) j 



P 



l+Tl 



(A.13) 



We present the coefficients of /z([3) in Tab. |3j For our choice of squark and gluino masses, m q 
400 GeV and mg = 500 GeV, the fit is accurate at the per mille level. However, the squark mass 
dependence of the coefficients at is weak, so that for the ratio rrig = 1.25m§, the fit may even be 
used in the whole range m q = 200GeV. . . ITeV with a an error of less than 1%, the maximum 
error being 0.7% at the lowest values of m q . 







,(10) 
J 88 




,(io) 

Jqq 




,(10) 


a\ 


A 

u. 


AC\A A A 1 11 
4911)1)1 12 


A 

u. 






,1442200/ 


a 2 


1 

— 1, 


A1 A A ^ £L A A 

91412644 


1 

7, 


76393459 


1 


63082127 


«3 


— (J 


6U19284U 


U, 


67354296 


A 
— (J 


23767546 


a A 


1 

1 


CCTi 1 Tin 

55221239 


—4, 


,45797953 


— (J 


59322814 


a 5 


-0 


00562902 





01846674 


-0, 


00167924 


a-6 





05091553 


-0, 


,16690138 


0, 


,01653721 


a-j 





,96279970 


-4, 


26865285 


-0, 


,88130118 


«8 


-0 


,32470009 





61779585 





07212666 


ag 


-0, 


12126013 


1. 


09469670 





,36722998 


aio 


-0, 


,50416435 





,85492550 


0, 


03809051 


a n 





70967285 


-0, 


,61743437 


-0, 


48810751 


a n 


3. 


99890940 


-8, 


,15404313 


-0, 


,02789299 


A13 


-9, 


83101731 


4, 


,99311643 


2, 


37790285 


«14 


-16, 


88306230 


44, 


33618984 


-0, 


60307995 


a\5 


52, 


83111706 


-26, 


25150811 


-7, 


,30816582 


«16 


20, 


,74447249 


-93, 


67735527 


1. 


.80344717 


an 


-99, 


69494837 


57, 


,74025069 


9 


30489420 



Table 3: Coefficients for the fit functions ©, QSJ and (lAlHl for m§ = 400GeV and m f = 500GeV. 

Note that the values for the constant a \ for quark-anti-quark scattering in Tab.[3]can also be used 
for other ratios of r q = mg/m q to very good accuracy because of the weak squark mass dependence 
at NLO. Since a\ = a\ (r q ), we can be extract from Tab.[3]for the value of a\ at r q = 1.25: 



< x _ 1681^_^__ 
ai[rq) ~ 400 (l+r2) 2 



a\(r q = 1.25). 



(A. 14) 



This is due to limp^ /?(|3)/|3 = a\ and Eq. (120 
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